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A different lattice fermion metliod is introduced. Staggered domain wall fermions are defined 
in 2n + f dimensions and describe 2" flavors of light lattice fermions with exact U(l)xU(f) chiral 
symmetry in 2n dimensions. As the size of the extra dimension becomes large, 2" chiral flavors with 
the same chiral charge are expected to be localized on each boundary and the full SU(2")xSU(2") 
flavor chiral symmetry is expected to be recovered. SDWF give a different perspective into the 
inherent flavor mixing of lattice fermions and by design present an advantage for numerical simu- 
lations of lattice QCD thermodynamics. The chiral and topological index properties of the SDWF 
Dirac operator are investigated. And, there is a surprise ending... 
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I. INTRODUCTION 

Lattice fermions are elusive. They not only present enormous challenges to numerical simulations of lattice QCD 
and other strongly interacting field theories but also pose in a most blatant way the problem of non-perturbative 
regularization of chiral gauge theories. Obviously the two problems have a common source. In the past several years 
enormous progress has been made in this direction. Interestingly, extra dimensions have been used again in theoretical 
physics, except this time the dimension is a tool to generate the correct low energy physics. 



Domain wall fermions were introduced in 0, | 



B 1^ 5| A large volume of work has followed since. The reader is 

referred to the annual reviews 

I' I' liiP'Tl @ references therein. The DWF latt ice regulator 

begins by defining a massive Wilson fermion ||17|| in 2n + 1 dimensions. If the boundary condition at the edges of the 
extra dimension is free then chiral surface states develop with the plus chirality fermion exponentially bound on one 
wall and the minus chirality fermion on the other wall The two chiralities have an overlap that breaks chiral 

symmetry. As the size Ls of the extra dimension increases the overlap tends to zero exponentially fast. As oo 
the theory has a single massless Dirac fermion in 2n dimensions. Obviously, this construction addresses both problems 
mentioned above. 

Since the DWF construction starts with massive Wilson fermions, it is easy to see that for any finite Ls there can be 
none of the exact chiral symmetries available in the naive lattice fermion formulation. In return for this shortcoming, 
flavor mixing between doubler states inherent in naive fermions are pushed up to the scale of the lattice cut-off, making 
them irrelevant. Early in the development of lattice field theories, the staggered approach to lattice fermions was 
devised to preserve some of the exact chiral symmetries of naive fermions. This meant there was still flavor mixing 
between the remaining doubler states |2^. In the spirit of staggered fermions, our domain wall construction 

presented here will preserve some exact chiral symmetry at any at the cost of introducing flavor mixing between 
light fermion doublers. Flavor symmetry violations should be exponentially suppressed in the — > oo limit. A 
preliminary version of this work was presented in [ p2[ . 

Staggered domain wall fermions (SDWF) are similar to DWF in the use of exponential localization of surface states 
to counter unwanted features of lattice regularization. In our case, this means disentangling the inherent flavor mixing 
between light doubler states. As a result, even at finite Ls, they have an exact U(l)xU(l) chiral symmetry very much 
like standard staggered fermions. This property makes them attractive for QCD thermodynamic simulations. As 
Ls ^ oo the light surface states in the theory are expected to recover the full SU(2")xSU(2") chiral symmetry. It 
must be noted that SDWF (and Wilson DWF as well) may not be able to develop light states if the coupling is 
extremely strong. Depending on one's perspective, SDWF combine the nice properties of the domain wall method 
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and staggered fermions. For simulations of QCD thermodynamics with standard DWF the reader is referred to 

We would like to draw the attention of the reader to a more subtle issue in this paper that might otherwise be 
overlooked. In most of what follows, the Saclay basis proposed by Kluberg-Stern et al. and predecessors 
is used for its nice spin-flavor algebra. While this is formally equivalent to the standard bases used for numerical 
simulations the transformation is gauge dependent and quite complicated. However, the conclusions we draw 

from our analytic work should be basis independent. The construction of actions suitable for numerical simulation 
will be discussed at the end of the paper. 

One interesting result is that it may be possible to do numerical simulations directly in the Saclay basis when 
Pauli-Villars fields are introduced. This applies equally to staggered fermions as well as SDWF and may eliminate a 
serious obstacle for using the Saclay basis for simulation. The reason is that the gauge field dependent transformation 
that accompanies the basis change will cancel as part of the subtraction. We believe this is yet another example of 
the potential usefulness of doubly regularized lattice fermions [ p5[ . 

For a semi-infinite extent in the extra dimension, the theory has four chiral fermions with the same chiral charges 
and is anomalous. To construct an anomaly free theory, we must use such "quadruplets" with charges as dictated 
by the corresponding anomaly cancellation condition. This is completely analogous to the case of Wilson DWF. Of 
course, in order to simulate a two flavor theory, a dynamical algorithm which effectively takes the square root of the 
fermionic determinant should be used. Further in the future, non-degenerate quark mass matrices could be explored 
to simulate the four lightest quarks. 

The paper is organized as follows: The SDWF Dirac operator and action is defined in section The symmetries 
associated with the SDWF action are presented in section [IF The flavor content is discussed in section |^ and 
the free propagator is calculated in secti on [v| The transfer matrix along the extra direction is given in section [v| . 
The promised surprise is in section VII but the reader ought to work through the preceding sections flrst... The 



transcription to the single component basis, problems and future directions are presented in section VIII, A discussion 
about alternative actions is given in section IX . The paper is concluded with section 



II. STAGGERED DOMAIN WALL FERMIONS 

In this section the SDWF Dirac operator and action is presented in the Saclay basis Here, we show that in 
the free theory light fermion fields localize exponentially along the extra direction with suppressed flavor mixing. 
The SDWF partition function is 

Z = J [dU] J [dWd^] J [d^U^] e-^. (1) 

Ufj,{x) is the gauge field, x is a site coordinate vector in the 2n dimensional space and jj, — 1,2, ■ ■ ■ , 2n. VE'(?/, s) is the 
fermion field and s) is a bosonic Pauli-Villars (PV) field, y is a hypercube coordinate vector related to the site 
vector xhyx = 2y + + A where O is a 2n dimensional binary vector indicating one of the 2^" possible origins of the 
hypercubic structure and A is a binary vector which indicates position within the given hypercube y. This implies 
the relations y = {x — O — A)/2 and A= [{x — O) mod 2] between the vectors, s = 0, 1, • • • , — 1 is a site coordinate 
in the 2n + 1 direction, where Ls is the number of sites in this dimension. 
The action S is given by 

5 = 5(/3,L„mo,m/;[/,'^,*,$)= (2) 
Sg{P; U) + SFiL„mo,mf,^, V) + Spy{Ls,mo; $^ $, V) 

where 

1 



-Re TrC/, 



(3) 



is the standard plaquette action with /3 — 2iVc/g§ with go the lattice gauge coupling and Nc the number of colors. 
The fermion action is 

Sp = - ^{y,s)DF{y,s;y',s')^{y',s') (4) 

with the fermion matrix given by 

Dpiy, s; y', s') = <5(s - s')D{y, y') + D^{s, s')S{y - y') (5) 



3 



where D{y,y') is the standard staggered action in the Saclay basis with the typical staggered mass (distance zero) 
set to zero and a different mass (distance one) proportional to (l/as — TOq) added as described below. Here are the 
expressions in a chiral basis in 2n = 4 dimensions. Extensions to other even dimensions are straightforward. 



D = 



B C 

-C^ -B 



(6) 



B 



mo 



^(10^5^) Am(I^) 



1 

2a5 



(7) 



(8) 



A^iV; y, y') = - [5{y + fi - y')V^{y) + 5{y-fi- y')V^{y') - 25{y ~ y')] , 



V^iV; y, y') = - [5{y + A - y')VM - 5{y - /I - y'W^iy')] 



(9) 
(10) 



where cri,2.3 are the Pauli matrices and is the identity. V^{y) are the gauge links between hypercubes related to the 
Ufi{x) links of the gauge action by Vfi{y) — U^{2y + 0) U^{2y + + jl). The parameter mo is the mass representing the 
"height" of the domain wall. The s-dependent part of the Dirac operator is exactly as for DWF but with a different 
mass mixing. Here we consider the action for 2" degenerate flavors. For non-degenerate flavors the action is similar 
and can be constructed according to the rules outlined in section [H . 



= D5 + Himf), 



(11) 



D5is,s') 



^ PRSil -s') s = 

i[PRS{s + l-s') + Pl5{s - 1 - s')] 0<s<Ls-l 
^ PlS{Ls -2-s') s = L,-l 



(12) 



The mass mixing term depends to whether Lg is even or odd. For odd Lg we have the following purely imaginary 
terms 



H{mf] s, s') 



1 

05 



imf [Pr5{s -Ls + 1)6{s') + PlS{s)5{Ls - 1 - s')] 



(13) 



where m/ is the degenerate mass of the flavor states localized on the domain wall. For even Ls more care must be 
used in constructing the mass mixing term to preserve the exact U(l)xU(l) symmetry (see section HI). In that case 
the mass term does not just involve the boundary at s = and L^ — l but also at s = 1 and — 2. Furthermore notice 
that the terms are real. The different "reality" of the mass term for Lg even/odd is just a reflection of a staggered 
wavefunction phase of the form The even Lg mass term is 



H{mf] s, s) — 



m/ 
as 



The chiral and flavor projectors are 

1±75 

Pr,L = P± ^ — - — , 



[PLBS{s)S{Lg -2-s')+ PllS{s - l)S{Lg - 1 - s') 
+PrrS{s - Lg + 2)S{s') + PrlS{s -Lg~ \)5{\ - s') 



Fr,l = F± = Prl,rl = P±± = (P± ® F±). 



(14) 



(15) 



The gamma matrices are taken in the standard chiral basis. For example 2ri = 4 dimensions they are chosen to be 

a 



7m 













74 = 




I) 


, 75 = 





(16) 
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with (Tp the Pauh matrices. The flavor matrices are defined as usual 



T 



(17) 



and notations like (75 ® ^5;^) mean (75 eg) ^zS,fj,). 

As with DWF, the PV action is designed to cancel the contribution of the heavy fermions. Thisjsjiecessary because 

531 . There is 



the number of heavy fermions is ~ Ls and in the Ls — > 00 limit they produce bulk type infinities ^, ^ | 
some flexibility in the definition of the PV action since different actions could have the same — > 00 limit. However, 
the choice of the PV action may affect the approach to the Lg 00 limit. Here the same approach as in |34, 3^ was 
chosen. The m/ = 1 case is exactly the quenched theory (infinitely massive fermions) . The PV action is 



(18) 



The symmetries and detailed properties of the SDWF Dirac operator will be discussed in the rest of the paper. 
However, as a first check we verify that in the free case the SDWF Dirac operator indeed describes four flavors with the 
chiralities localized on the opposite walls. Following identical steps as in ^ we go to momentum space and demands 
that in order for light modes to exist there must be a wavefunction such that 



Dpik, s; k', s')4>{k',s') = Aiaivc(fc, s- k', s')c^{k', s'). 



(19) 



In essence this equation demands that the extra term in the Dirac operator cancels the flavor breaking term B. 
From the above equations it is easy to see that Eq. ([l9| ) leads to 



^ [P+5{s + 1 - ,s') + P^5{s - 1 - s')] + E^(75 ® - s')] 0(fc, s') = 



(20) 



where 



1 — cos k^ — mo + 1/05 



(21) 



From this equation, it is easy to sec that the P± projectors in the s-dependent part commute with the flavor breaking 
part so that each may be simultaneously diagonalized. This constraint alone effectively restricts the allowed projectors 
to the ones chosen here. 

The solution is separable and 0(s) is the s-dependent part 



where 



where 



-+{s) — P ^4'{s), etc. In this notation, we can write 

/ 

V 



-) 



(22) 



-b 



(23) 



J = 1,2,3, 64 



(24) 



Solving Eqs. ( |20| ) relating nearest neighbor s sites is a bit complicated because the flavor components mix and is not 
of interest for this discussion. On the other hand the solutions to these equations after iterating twice are simple. For 
as = 1 we have 



4±+(s±2) = -66V±+(s), 
6±_(.s±2) = -Pb^i^is). 



(25) 



For free fermions, [6,6'] = and 66', 6'6 are both proportional to the identity with eigenvalue 



A(66^) = A(6^6) = 6^, 



(26) 
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'E,bl. (27) 
If we require that 

< 1 (28) 

then for a semi-infinite s direction, s>0, only (t>+± is normahzable, while 0-± is not. However, this is not enough to 
ensure that the doubler modes are not present. We must further require that the above condition excludes momenta 
with components larger or equal to tt. This can be seen by writing out Eq. ( p8|) 

= J^[(l-cosfc^) + (l-mo)]'<l. (29) 

For momenta near the origins of the n'^ Brillouin zone (where n is the number of momentum components near tt) 
this gives 



When n = the sufficient condition is 



(1 - mo)^ + n(l - mo) + n < 1. (30) 



< mo < 2, (31) 



same as for Wilson DWF. However, this condition does not ensure that the doubler modes are non-normalizable for 
all n. For example, for n — 1 and mo = 1.5 both Eqs. (|3^) and (|l]) are satisfied making the 1^* Brillouin zone doubler 
wave functions normalizable. The range of mo needs to be further restricted. The following condition ensures that 
only the 0'^ Brillouin zone wave function is normalizable 

< mo < 1. (32) 

The above is presented graphically in Fig. |l|. It is straightforward to extend these results to the more general case of 
< 05 < 1. 

III. SYMMETRIES 

When constructing the SDWF action, it is important to preserve the symmetries of the massless staggered action 
p9| , |3^ . Of course, adding any new terms to the staggered action will break some of those symmetries, so we have 
to find new symmetries that involve the extra dimension. The symmetry transformations for the action in the Saclay 
basis of section |l| are presented below. 

U(l)eXU(l)o chiral rotations. The presence of this symmetry is one of the main motivations of this paper. The 
residual chiral symmetry of staggered fermions involves making separate chiral rotations on even and odd sites. 
Terms in D-^ are not invariant under these rotations unless we extend the notion of even and odd, including the extra 
dimension. The operator iS is defined as 

Ss^s' = i-iySis - s') (33) 
and then the extended even/odd projection operators are defined as 

Pe = i[(l®l)+ 5(75 0^5)], 

Po = i[(l®l)- 5(75 0^5)]. (34) 

Using these projection operators the chiral transformation is 

*(?/) ^ (e^'^Pe + e^'^Po) *(?/), 

*(2/) ^ *(y) (e^'^Pe + e^'-Po) ■ (35) 



Rotations by n/2. These rotations are in planes perpendicular to the extra dimension and the transformations are 
the same as the original staggered ones. 



a=l, Z (l-cos/7 )/2 = n 



1 1 r 




0.0, 



FIG. 1: From Eq. (§9): b'^ at the origins of the five Brillouin zones n — [0, 1, 2, 3, 4] plotted vs. mo. 



fi-parity. These transformations reflect the 2n — 1 spacetime axes perpendicular to the spacetime axis in the fi 
direction. is not invariant under this symmetry unless we also reflect the s direction as well. If the reflection 
operator is defined as 



then the transformation is 



7^. 



s) 

*(y,s) 



6{L, 



1 



^'(y,s') 7^,^,(7^®C5)• 



(36) 



(37) 



Shift by one lattice spacing. The toq — l/as term in the SDWF action breaks this standard staggered fermion 
symmetry at the expense of absorbing the renormalization of the flavor breaking term. That these terms are additively 
renormalized in the interacting theory follows from the work of Mitr a and Weisz . Nevertheless, interesting methods 
to alleviate the breaking of this symmetry are discussed in section VII] . This symmetry relates to interactions inside 
a hypercube which are essentially non-physical. We feel that the breaking of this symmetry is a small sacrifice, but 
the issue certainly can and should be debated. 

The symmetry transformation for the case toq = 1/^5 is given. Already for staggered fermions, the symmetry 
transformation is complicated in the Saclay basis due to the imposed hypercubic structure of the formulation. For 
SDWF, there is an added complication. Some parts of the transformation require a reflection in the s direction 

^{y) i [(1 - (7^5 ® i5)n] ^{y) 



+i [(1 ® Cp) + (7m5 C5)7^] «'(y + /t). 



*(2/)-*(2/)H(l®Cp) 
+'^(2/ + A)H(l®00 



^(75m ®^5)] 

-7e(75^®e5)] 



(38) 
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where s indices have been suppressed. 



IV. FLAVORS OF SDWF 



In this section the SDWF flavor identification is made in the Saclay basis. From section || we see that for a finite 
extra direction with Lg sites the components of aU flavors are locahzed around s—0 while the P_ components 
are locahzed around s=Ls — l. However, as already mentioned in section III one of the main goals of this paper 
is to preserve most of the staggered symmetries and particularly the U(l)eXU(l)o chiral symmetry. For example, 
to generate the four-dimensional flavor components q with we should choose s near zero. If s=0 is chosen, 

P++q{y) = P++\I'(j/, 0), then these components also belong to the Pe part of the fermion fleld. Therefore, to project 

flavor components with P ^ we are not only restricted to choose s near L^ — 1 but also choose s so these components 

belong to the Po part of the fermion field. Then, components P±+q will not mix even for finite Lg because of the 

even/odd symmetry. In this example, we would like to pick P i-q{y) = P y'^iy, s) with s being even and near Lg — l. 

So, if Lg is odd then s=Ls — l is a good choice. However, if Lg is even, then we should choose s^Lj— 2 instead. 

/X 

X 
X 
\ X , 

fx 
X 

X 
X 



( X 

X 

X 

\x , 

(X 

X 
X 

Vx 





(39) 



Using the block notation of Eq. (|22|), an example for even Lg is sketched in Eq. (|39|). In this equation ^'(s=0) is at 
the top and '^{s=Lg—\) is at the bottom. The capital letters denote one of the correct choices. On the other hand if 
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Lg is odd, e.g. Ls=S, then 



(40) 



We note that other choices for selecting flavor components near the boundaries are certainly possible. 



THE SDWF PROPAGATOR 



The propagator in the Saclay basis in momentum space and for my has the form 

D-\s, s') = Gi e(s - s') + Gs e{s - s' - 1). 
For m-f 7^ and Ls — multiple of 4 has the form 

D-\.s, s') = [Gi + m/Ga] e(s - s') + G3 e(s - s' - 1). 



(41) 



(42) 



For Ls — even but not multiple of 4, the propagator has the same form as in Eq. (^) but with ruf ^ ~^f- Fo^' 



Ls — odd the mj term is more complicated and is not given here. Nevertheless, based on sections |I| and III one can 
see that the form is similar. 

Giip; s, s'), G2(p; s, s') and G3(p; s, s') (p is the momentum) are proportional to the identity in their flavor indices. 
Also, Gi anti-commutes with (75 ® 1) while G2 and G3 commute with (75 (g) 1). The flavor mixing is in the function 

e(s) 



e{s) = (1 (?) 1), (s even) 



-, (s odd) 



(43) 



where |6| is given in Eq. (p7|). For s — s' even and 'mf=0 the propagator has no flavor mixing except for the G3 
term in Eq. (^l]). In this term e(odd) breaks flavor in exactly the same way as for free staggered fermions. An exact 
U(l)xU(l) symmetry is maintained. The matrix coefficient G3 vanishes exponentially fast with Ls for s, s' near 
opposing boundaries and therefore as Ls ^ 00 with m/=0 the propagator anti-commutes with (75 (E) 1) and has no 
flavor mixing provided s — s' is even. This is in accordance with the discussion in section For s — s' odd more 
severe flavor mixing is present. 
The Gi, G2 and G3 terms are 



Gi(p : s, s') = ;^7^, sinp^ [P+G(s, s') + P-G{Ls - 1 - s, P, - 1 - s')] 



(44) 



G2 = [P+Sis, 0)G{Ls ~ 1, s') ~ 5(8, Ls - l)P-G{Ls - 1, P^ - 1 - s')] 



(45) 



G3 = i^'-''-^^P+[hG{s,s') + ie+{s)G{.s-l,s')] 

+ P_ [hG(Ls - 1 - s, Ps - 1 - s') + ie-{s)G{Ls - 2 - s, - 1 - s')] 



(46) 
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where 

0+is) = 1, (s^O) 

0+{s) = 0, (,s = 0) 

0^{s) = 1, {s^L,-l) ^^'> 

0-{s) = 0, {s = Ls-l). 

There is an exact correspondence with the terms of the standard DWF propagator. In the notation of ^ the 
symbol correspondences between SDWF and DWF are 

G{p:s,s')^G+{p:s,s'), b ^ b, P^\p. (48) 

The reader is referred there for the detailed forms. 

As can be seen the decay coefficient is now in terms of b given in Eq. ( p7| ) instead of the DWF b = ^^[^ — cosp^] + 
l/a5 — mo. The localization condition for toq is over half the DWF range (see Eq. (|3l|)). Also, for Ls odd, the effective 
mass TOeff has the same general form as in Wilson DWF ^ 

mcff = (1 - Y™o)(^"/ + |1 -"io|^=). (49) 



VI. THE SDWF TRANSFER MATRIX 



The SDWF transfer matrix in the Saclay basis is presented. We can use the technique of Neuberger to rewrite 
the free SDWF determinant in a form that allows for a quick identification of the transfer matrix. A complete 
Hamiltonian analysis is beyond the scope of this work. 

After interchanging various rows and columns of the SDWF matrix, the determinant is equivalent to the determinant 
of the matrix 



Pi 

V 



/5o \ 



(50) 



where all of the as and j3s are the block triangular matrices 



-B 
C l/as 



, 13s 



B 



(51) 



For aL^-i and /3o, l/os is replaced with — /x/os so /i is a parameter that controls the boundary conditions: ^ = ±1 
for (anti)periodic and /i = for free. This parameter is used only for the derivation of the transfer matrix and it 
is not needed in the theory. In particular, the reader is cautioned against using /i as a mass parameter since it is 
inconsistent with the rules of section IV. The definitions of B,C are given in Eqs. (0) and (|^). 



In this notation, following Neuberger's construction leads to the SDWF determinant 



Df = (fixed sign) (det 5/05)^-' det 
and the SDWF transfer matrix identification 

T = 

We can easily check that 



-/i \ _ L, / 1 

1 j V -/^ 



5-1/05 B-^C 
C^B-^ ar^[C^B-^C - B] 



[^,(10^5)] = 0, {5,(10^5)1 = 



(52) 



(53) 



(54) 



and as a result 



{r,(i®e5)} = o. 



(55) 
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Also since = we can see that T is also anti-Hermitian 

= -T. (56) 

This is different from Wilson DWF and gives some idea why solving the zero mode problem in Eq. (|25| ) simplifies 
when solving for the field two sites away. Obviously, standard transfer matrix manipulations should be done with the 
Hermitian transfer matrix which corresponds to a Hermitian Hamiltonian H. 



VII. ABOUT THAT SURPRISE... 



At this point the reader must be wondering: "What is the spectrum of the transfer matrix?" and "What is the 
corresponding Hamiltonian?" Here is where we were a bit surprised. Just in case the reader will later (after reading 
this section) be tempted to claim that there is no surprise, she/he is invited to guess the 05 — > oo limit Hamiltonian 
as well as the general 05 spectrum of T^. 

We find, after some algebra, the following Hamiltonian Hq corresponding to the 05 — > limit of the transfer matrix 

. i E„ + ...0 ^ c ^ ^ 

The first surprise is that Hq is exactly diagonal in fiavor. It is almost, but not exactly, the same as the standard 
overlap Hamiltonian Hu,, as has a factor of ^ compared to C. Nevertheless, because Hq is diagonal in fiavor, the 
standard machinery of DWF can be directly applied. 

For example, let us consider the spectrum of Hq- Following methods identical to ||, ^ ^ we find that 

Ho(^^^=0 utAu + wtAw - 2mo 0. (58) 

Because the matrices S{y + /t — y')Vp,{y) are unitary the range for which this equation can have a solution (for all 
Brillouin zones) is 

< mo < 2 (as ^ 0). (59) 

For mo in the above range Hq can have zero eigenvalues that via the overlap formalism correspond to a change of 
index and to exact and robust zeros of the fermionic determinant. This can also be seen graphically in Fig. ^. The 
background field configuration is a smooth gauge field that has non-trivial topology |3^ . The plaquette value for 
that configuration (i.e. the sum in Eq. (^)) is about 0.05. For the rest of the paper we will refer to this configuration 
as the "instanton" configuration. 

The crossing diagrams were done for fermions in the Saclay basis with gauge fields VJ^ (y) defined on a 2^ lattice of 
hypercubes. Since no attempt is made here to make a connection with the topological charge of the underlying gauge 
configurations U^{x), the reader can regard the gauge fields as examples "pulled out of a hat." All numerical analysis 
was done by full diagonalization of the relevant matrices using the LAPACK libraries and an IBM-T20 Think Pad 
(which by the way performed brilliantly). 

However, numerical simulations are done at non-zero 05, typically at 05 = 1. In standard DWF there is an exact 
connection between zero eigenvalues of Hq and unit eigenvalues of the transfer matrix at any 05. This correspondence 
does not hold here. As a result the analysis is more complicated. In other words 

|eig(-r2)| = l ^ cig(i7o) = 0, as^O. (60) 



From Eqs. ( p5[ ) and (56) we deduce that the spectrum of — T is strictly real, positive and doubly degenerate because 
T is anti-Hermitian and anticommutes with (1 0^5). In 2n = 2 dimensions that would be the end of the story because 
there are only 2" — 2 flavors. Two exact zero modes are produced for every crossing in the Hamiltonian spectrum. In 
2n = 4 dimensions there are 2" — 4 flavors and we do not have an exact correspondence between the degeneracy of 
the unit magnitude eigenvalues of T and the number of flavors. This is problematic since this is a basic and defining 
property for a chiral theory. 

For the "instanton" background the eigenvalue crossing diagram of 

\tm = log[spectrum (-T^)] (61) 

vs. tuq is given in Fig. ^ for 05 = 1. A closeup of the region around mo=0 is shown and the eigenvalues are marked 
depending on the sign of the 9(75 ® !)(/ where q is the corresponding eigenvector. All eigenvalues in Fig. |^ are four-fold 
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m 



FIG. 2: The spectrum of Ho vs. mo for an SU(3) "instanton" background on a 2"^ lattice of hypercubes. 

degenerate! The symbols for each eigenvalue are there but are literally overlapping. An example is given in Table |. 
Also, the reader should observe that in Fig. || the transfer matrix eigenvector chirality does follow each flow line. This 
ensures that the four modes are of the same chirality and therefore a crossing should correspond to a net change of 
four in the index. It is of fundamental importance that this occurs here as it reassures us that on a given boundary 
of the extra dimension there are four flavors of light chiral fermions with the same chiral charge. 

For a very rough background gauge field (plaquette 0.85) the crossing diagram is given in Fig. ^. The eigenvalues 
are doubly degenerate and very nearly but not quite four-fold degenerate. The non degeneracy is small and almost 
non visible. One of the worst cases is presented in Table ||. The lack of exact four-fold degeneracy should be the 
subject of further research but it is obviously small even on this extreme background gauge field that is unlikely to 
occur in current numerical simulations of QCD. Furthermore, it would be interesting to study the transfer matrix in 
2n > 6 dimensions, where the required near-degeneracy should be 2", to see if the same behavior persists. 

Even with the lack of exact four-fold degeneracy, if we choose mp away from the crossing region then the fermion 
determinant will still have four exact zeros. This will only break down if mp is chosen between the two nearly 
degenerate sets of crossings. Of course, based on DWF studies at lattice spacings used in today's simulations one 



12 



0.4 



0.2 



0.0 



-0.2 



-0.4 




0.0 



0.2 



0,4 



m 







FIG. 3: Xtm vs. mo for 05 = 1 and an SU(3) "instanton" background (plaquette ~ 0.05) on a 2"^ lattice of hypercubes. The 
diamonds represent eigenvectors of with chirality +1 while the pluses with —1. All eigenvalues are four- fold degenerate and 
are indistinguishable by the graphics. For an example please refer to Table pi 



TABLE I: The near zero spectrum of log(— T^) for mo = 0.2 for an SU(3) "instanton" background on a 2"^ lattice of hypercubes. 



mo 


logA(-T2) 


0.2 


-0.0408683 


0.2 


-0.0408683 


0.2 


-0.0408683 


0.2 


-0.0408683 


0.2 


0.327367 


0.2 


0.327367 


0.2 


0.327367 


0.2 


0.327367 
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FIG. 4: Aim vs. mo for — 1 and a "rough" SU(3) background (plaquette ~ 0.85) on a 2* lattice of hypercubes. All eigenvalues 
are two-fold degenerate and are indistinguishable by the graphics. Although there is no exact four-fold degeneracy the graphics 
can hardly distinguish the non-degeneracy. Please refer to Table ^ 



expects dense crossings in the usable range of mg. Then there will be configurations for which uiq is between double- 
crossings that have split. As far as topology is concerned this will break flavor to some degree. Nevertheless, since 
in DWF the dense crossings correspond to small instantons and are unphysical we would expect that as the lattice 
spacing becomes smaller such configurations will become less important. 

Nevertheless it is instructive to further study the spectrum of T for unit magnitude eigenvalue. Because T contains 
the inverse of B it is hard to study analytically. However, in the subspace \Xtm\ = 1 we can proceed as in [||, ^, 
In particular we find that 

t( =±tT\ - \ => Hp \ - \=0 (62) 
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TABLE II: The near zero spectrum of log(— T^) for mo = 0.3 for a "rough" SU(3) background (plaquette ~ 0.85) on a 2'* lattice 
of hypercubes. 



TtlQ 


Luy/\y 1 ) 


u.o 


u.^yuu^iij 


u.o 


n ^Qond^ 

u.^yuu^o 


u.o 




u.o 


n 9Sfi9Ql 


u.o 


u.^u^uoo 


u.o 


u.^u^uoo 


u.o 


n 9fil 'i'll 


u.o 


n 9fii 'I'li 


u.o 


n n9fin'ifiR 

U.U^UUOUU 


0.3 




0.3 


0.0272779 


0.3 


0.0272779 


0.3 


0.414225 


0.3 


0.414225 


0.3 


0.414597 


0.3 


0.414597 



where 

So, we can study the crossings of the spectrum of this pseudo-Hamiltonian Hp. The crossing range can be determined 
as before by using the unitarity of the matrices 5{y + p- — y')V^j,{y). There are two crossing ranges 




2 2 

< mo < 2 and — < 7710 < h 2. (64) 

as as 



For the "instanton" background the crossing diagrams (for aU BriUouin zones) are shown for 05 = 1 in Fig. ^ 05 = 0.5 
in Fig. ^ and 05 = 0.25 in Fig. ^ and the reader can see the agreement with Eq. (|6^). 

But this is not aU. We had a hard time at first because we used SU(2) gauge fields. And the spectrum was always 
inexplicably not two- fold but four-fold degenerate, for any gauge field (smooth or rough)! This is the second surprise. 
And it will not be investigated further here. After all, what good is a paper that does not leave some mystery behind... 
In any case this property is likely linked to the fact that —1 is part of the groups SU(2), SU(4), . . . but is not part of 
the group SU(3), etc. 



VIII. THE SINGLE COMPONENT BASIS AND SIMULATING SDWF 



In the previous sections we discussed the properties of the SDWF Dirac operator in the basis proposed by the 
Saclay group. By now, the reader has seen the advantage of this approach in determining the properties of robust 
zeros of the transfer matrix. In the past, technical problems have disfavored direct simulation of dynamical fermions 
in this basis in lieu of the simpler single component basis, where there is just a single fermionic spin-flavor degree of 
freedom per site. 

In this section, we will first discuss the construction of the SDWF in the single com ponent basis with an emphasis 



on preserving, to the maximum extent possible, all of the symmetries of section III. One consequence is that the 
spin-flavor algebra, i.e. (75 (g) 1)^ = (1 (g) 1), will be broken and only recovered in the continuum limit. Next, we will 
discuss a technique where the spin-flavor algebra is restored at the expense of some lattice symmetries. We believe 
these lattice symmetries can be restored by proper stochastic averaging. Finally, we will discuss a new algorithm 
based on double regularization for simulating staggered (and SDWF) fermions directly in the Saclay basis. As is 
typically the case in our fleld, performance during numerical simulation of QCD will likely determine which of the 
three proposals survive. We feel that further research in this area is needed. 
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FIG. 5: Xtm vs. mo for as — 1 and an SU(3) "instanton" background (plaquette ~ 0.05) on a 2"^ lattice of hypercubes. 



A. Single component basis without projection 



Let us first review the staggered action in the single component basis and the transformation connecting it to the 
Saclay basis. From Eqs. (|^) through (p^, with mp = l/as, the equivalent 2n dimensional Dirac operator is usually 
written 



D{x, x')x{x') = \ ^(-l)A-')(-) [U,{x)5{x + A - x') - Ul{x)5{x - A - x')] x(x') 



(65) 



with the components of the binary vector rj given by 



?7i(x) = 0, 

772(2:) = xi (mod 2), 
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FIG. 6: Xtm vs. mo for as — 0.5 and an SU(3) "instanton" background (plaquette « 0.05) on a 2* lattice of hypercubes. 



mn{x) = xi -\ ^ X2n-i (niod2). 



(66) 



In the free theory, the unitary transformation from single component site-wise fields x(x) to the hypercubic fields 
ilj{y) is simple. As in section ||, if we label sites on the hypercube starting at the origin O by a binary vector A the 
transformation is 



^ aa, 

Axi^y + o + A) 



(67) 



with the rows of the 2^" x 2^" dimensional matrix T indexed by the various combinations of spin a and flavor a indices 
and the columns indexed by the corners of the hypercube y. Specifically, the components of the free F may be chosen 
as 



(68) 



where the row index a and column index a of the 2" dimensional representation of the Clifford algebra are interpreted 
as spin and flavor indices. 
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FIG. 7: Xtr, 



mo for as = 0.25 and an SU(3) "instanton" background (plaquette « 0.05) on a 2* lattice of hypercubes. 



When the staggered action is made gauge invariant, differences will arise between the two formulations. In the 
single component basis, fermions are only coupled to nearest neighbor sites, so simple links are all that is needed to 
preserve gauge invariance. Of course, longer paths could be used to link sites and indeed are often used to implement 
an improvement program. In the Saclay basis, since fermion fields are associated with hypercubes, the gauge fields 
must be used to move the components on the hypercube to some common point where the hypercubic field can be 
assembled. While it may seem that the resulting actions could be completely different, they have the same terms, 
site by site, that differ only in the choice of paths used to connect nearest neighbors. Hence, they have the same 
continuum limit. 

However, this is not the end of the story. For example, when moving a component two sites on the hypercube for the 
construction of the hypercubic field there are two equivalent minimum distance paths from which to choose. Choosing 
one path over the other will preserve the unitarity of F but break the rotation by 7r/2 symmetry. Conversely, choosing 
to average over both paths preserves rotations but means F need not be unitary, and thus potentially singular for 
sufficiently rough gauge fields. For SDWF, some terms in the action may also break the shift by one lattice spacing 
symmetry. The well known source of the problem is the imposition of the artificial hypercubic structure for the 
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identification of spin and flavor degrees of freedom ||3^ . 

Thus, a conservative approach is to abandon a transcription from hypercubic bases and directly use the single 

component basis and the technique of Golterman and Smit |^ . Using the symmetric shift operator with the binary 
vector ({x) and scalar e{x) 

E^{x,x') = (-l)'i<(-)i [U^{x)6{x + fi^x') + C/t(x - - A - x')] (69) 

Cl{x) =X2^ h X2n (mod 2), • • ■ , (2n-l = (mod 2), C2n = (70) 

e{x):^xi-\ \-X2n (mod 2) (71) 



we write the toq term 



^ - E ® ^5;.) V'(y) ^ - "^o) (-l)^(^) E ^m(^> ^')x(x') (72) 

and the chiral projection operators are constructed using 



(75®a)^(y)-^YSr ^ eMi-M.„^m(^,^''')x---xi?^,„(a:(2"-i),x(2"))x(x(2")) (73) 

where e/ji - /i2n totally antisymmetric tensor and the summation over the 2n site vectors a;*-^^ , • • • , x^^"-* is implied. 

The upside to this approach is that it preserves the staggered symmetries to the extent possible. The downside is 
that chiral projection is no longer exact, except in the continuum limit. Of course, this is a different manifestation of 
the same problem that makes the transformation to the Saclay basis non-unitary, where projection is exact. 



B. Exact projection with stochastic symmetrization 



The second approach addresses the projection problem at the expense of breaking some symmetries, which can be 
restored in the ensemble average as described below Our example will use the hypercubic basis of Daniel and 
Sheard | ]3l[ but equivalent examples are to construct a unitary transformation T into the Saclay basis or even to 
restrict the Golterman-Smit operators to single paths between sites. 

Quickly reviewing the Daniel-Sheard formulation, we want to construct "local" fermion bilinears of definite spin and 
flavor, where local means local to the hypercube, from the single component states xi^)- We identify the hypercubic 
Daniel-Sheard fields by a simple relabeling: XA{y) — xi^) ^-nd x — 2y + O + A as before. Local bilinears are written 

XAiy)hF®l^)ABXBiy) = Y.i-if^^-^'HxMx, x')x{x') (74) 

where x,j;' are summed over the hypercube and lA{x^x') represents the hnks chosen to make the bihnear gauge 

s s 

invariant. The notation is 75 = ^ x ■ ■ • x ^27^ phase factor is computed from 

{is <S> ^f)ab => HA B) = ^Tr (7^75757!^) • (75) 

As an aside, this gives exactly the same terms appearing in Eqs. (^2|) and ( |7^ ) provided you keep only the terms on 
a single hypercube. 

As mentioned before, imposing a hypercubic structure introduces problems with maintaining the staggered sym- 
metries for arbitrary spin and flavor choices. Our proposal is at the beginning of each update step of whatever 
update algorithm, first choose the origin O at random from the 2^" ways of imposing the hypercubic structure on 
the lattice. Next, choose at random only one of the minimum distance paths on the hypercube for making bilinears 
gauge invariant with the restriction that the same path is used in both directions. Thus, U{A,B) is unitary and 
W^iA, B) = U{B, A). Note that different paths can be used on different hypercubes. Choosing a random hypercubic 
structure and random paths on the hypercubes at each update step ensures that symmetry breaking effects due to 
these choices will cancel out in the ensemble average. The purpose of choosing only one path per pair of corners on 
the hypercube is to guarantee the chiral projection property. For example, (75 ® 1)^ [(75 ® 1)^] —f x which is the 
same as the continuum where we normally choose Hermitian gamma matrices. 
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C. Doubly regularized staggered fermions 



The third proposal is specific to the Saclay basis but apphes equaUy weU to SDWF and staggered fermions with 
Pauh-Villars fields. Following the second proposal, we can construct at each update step a unitary transformation 
from the single component basis to the Saclay basis that will certainly depend on the gauge field but not on the 
value of the mass m/. Since the fermionic action and the Pauli-Villars action only differ by the value of m/, then 
the contributions from the transformation will cancel between fermions and the pseudofermions. Specifically, the 
fermionic partition function on a fixed gauge background and for finite lattice spacing, volume and is 

gXFr^CF(m^)rxF-0^.vr+-DF(m^=i)r0pv 

_ detr'l' deti:>F(m/) detT _ detDpCm/) 
~ detrt c\elDY{mf ==1) detT ^ detDplm/ = 1) 

= j [d^d^ j [d^U<^] e*^F(™/)*-*'^F('"/=i)* (76) 

which is what we had back in Eq. (]^). In practice, we do not need to specify the paths chosen for the basis transfor- 
mation since they cancel from the path integral. But, it would still be important to choose at random the origin O 
of the hypercubic structure at each update step to avoid violations of the shift by one lattice spacing symmetry. We 
would like to emphasize again that this proposal should work for staggered fermions with added Pauli-Villars fields 
and we believe this is another example of the potential of double regularization to improve the usefulness of existing 
fermion actions by canceling lattice artifacts |2^. Also, notice that the number of fermionic degrees of freedom is the 
same as in the single component basis because the Saclay fields are defined on hypercubes. 

Of course, some of the ideas presented in this section for implementing SDWF for numerical simulation have been 
discussed before, e.g. the idea for stochastic restoration of staggered symmetries is certainly descended from the work 
of Christ, Freidberg and Lee p9| |. 
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[c^XfC^Xf] / d(j)l,yd<j)-_ 



ypv 



IX. ALTERNATIVE ACTIONS 



The SDWF action considered here is not unique. It is possible that actions with better scaling properties may be 
constructed using improved fields in the same spirit as with staggered fermions (see 4^ and references therein). 
Additionally, in our earlier work we introduced the domain wall defect using a local mass term (distance zero) 
which preserved the shift by one lattice spacing symmetry (among others) and broke the U(l)xU(l) chiral symmetry. 
In this work, we considered a distance one mass term which preserves the chiral symmetry and breaks the shift 
symmetry. We view this as a better choice because the additive renormalization it produces ]37| merely contributes to 
the flavor breaking term that our domain wall formulation is designed to eliminate. It is possible that other distance 
mass terms might prove useful in the future and even have a faster exponential rate of restoration of flavor symmetry. 
We emphasize that the primary requirements for these mass terms are that they be of the order of the cutoff and 
commute with the operators in Eq. (20). 



X. CONCLUSIONS 



In this paper a different lattice fermion regulator was presented. Staggered domain wall fermions are defined in 
2ji + 1 dimensions and describe 2" flavors of light lattice fermions with exact U(l)xU(l) chiral symmetry in 2n 
dimensions. The full SU(2") xSU(2") flavor symmetry is recovered as the size of the extra dimension is increased. 
SDWF give a different perspective into the inherent flavor mixing of lattice fermions and by design present an 
advantage for numerical simulations of lattice QCD thermodynamics. We have paid particular attention to the chiral 
and topological index properties of the SDWF Dirac operator and its associated transfer matrix. In the limit whore 
the lattice spacing in the extra dimension as tends to zero the corresponding Hamiltonian Hq is proportional to the 
identity in flavor space illustrating the complete absence of flavor mixing. 

For a semi-infinite extent in the extra dimension, the theory has four chiral fermions with the same chiral charges 
and is anomalous. To construct an anomaly free theory, we must use such "quadruplets" with charges as dictated by 
the corresponding anomaly cancellation condition. This is completely analogous to the case of Wilson DWF. 

However, there are still a number of unresolved issues related to this formulation which need to be studied in future 
work. In particular: 
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1 For QCD, the nearly four- fold crossing degeneracy of the Hamiltonian must be investigated thoroughly. 



2 SDWF should be implemented for numerical simulation according to the proposals of section |VIII[ . In the broken 
phase of QCD, it is obviously important to confirm that one pion is a pseudo-Goldstone boson and that the 
remaining fourteen non-singlet pions become degenerate with the pseudo-Goldstone boson as La —>■ oo. Also, 
the expected robustness of topological zero modes should be confirmed as was done for DWF |Q . 

3 We have presented an analysis of the zeros of the SDWF Hamiltonian through the pseudo-Hamiltonian Hp and 
in the limit 05 — > where the flavor breaking is trivially absent in Hq. A derivation and analysis of the full 
spectrum of the Hamiltonian for general as is needed. 

4 Since the nearly degenerate four-fold crossings in the spectrum of the Hamiltonian have the same chiral charge, 
the conserved currents of the full SU(2")xSU(2") symmetry must exist and can be constructed in the overlap 
formalism p^ . Simpler constructions of these currents may be possible. Constructing and measuring the 
conservation of these currents in simulations is important. 

5 SDWF were constructed with the simulation of QCD thermodynamics in mind because of the importance of 
having a continuous subgroup of chiral symmetry for any Ls . It is worth confirming that this gives SDWF some 
advantage over DWF in looking for critical fluctuations at the finite temperature QCD phase transition. 

6 It would be very interesting to add the Kogut-Sinclair four-fermion interaction Q to the SDWF action to 
enable simulation at zero quark mass in the region of the QCD phase transition. 

7 The domain wall mass term must be of the order of the lattice spacing and will introduce a hard breaking of 
some part of the staggered symmetry group, in our case the shift by one lattice spacing symmetry, causing 
quantum corrections to the SDWF transfer matrix |^ . Our analysis of the transfer matrix spectrum indicates 
a range of mo values can still be found where the transfer matrix behaves correctly in the presence of gauge 
field topology, as was the case with Wilson DWF. We believe this issue should be studied thoroughly ||4^. 

8 Can SDWF reveal (or has it already revealed) some new insight into the nature of inherent flavor mixing of 
lattice fermions? 

On the strength of the results of our transfer matrix analysis, we believe that SDWF may be an attractive alternative 
to DWF. The formulation is now sufficiently mature that the issues above should now be addressed in the context of 
QCD. As has been the case with staggered and Wilson fermions in the past, we should have a choice between SDWF 
and DWF for a given problem according to our resources and preference in the near future. 
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